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O 1 ABSTRACT 

We present the results of the first long-term (2.2 years) spectroscopic monitoring of a gravitationally lensed quasar, namely the Einstein Cross 
^ , QSO 2237+0305. The goal of this paper is to present the observational facts to be compared in follow-up papers with theoretical models to 
constrain the inner structure of the source quasar, 
(f-) . We spatially deconvolve deep VLT/FORS 1 spectra to accurately separate the spectrum of the lensing galaxy from the spectra of the quasar 
^ ' images. Accurate cross-calibration of the 58 observations at 31-epoch from October 2004 to December 2006 is carried out with non- variable 
GO foreground stars observed simultaneously with the quasar. The quasar spectra are further decomposed into a continuum component and several 
broad emission lines to infer the variations of these spectral components. 

We find prominent microlensing events in the quasar images A and B, while images C and D are almost quiescent on a timescale of a few 
months. The strongest variations are observed in the continuum of image A. Their amplitude is larger in the blue (0.7 mag) than in the red 
ON ' (0-5 mag), consistent with microlensing of an accretion disk. Variations in the intensity and profile of the broad emission lines are also reported, 
CD most prominently in the wings of the C III] and center of the C IV emission lines. During a strong microlensing episode observed in June 
t^- , 2006 in quasar image A, the broad component of the C III] is more highly magnified than the narrow component. In addition, the emission 
lines with higher ionization potentials are more magnified than the lines with lower ionization potentials, consistent with the results obtained 
with reverberation mapping. Finally, we find that the V-band differential extinction by the lens, between the quasar images, is in the range 
0.1 -0.3 mag. 
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1. Introduction 

The gravitational lens QSO 2237+0305, also known as 
"Huchra's lens" or the "Einstein Cross", was discovered by 
Huchra et al. d 1 985t during the Center for Astrophysics 
Redshift Survey. It consists of a z s = 1.695 quasar gravita- 
tionally lensed into four images arranged in a crosslike pattern 
around the nucleus of a zi = 0.0394 barred Sab galaxy. The av- 
erage projected distance of the images from the lens center is 
0.9". 

A few years after this discovery, Schneider et al. (1988) 
and Kent & Falco (1988) computed the first simple models of 
the system, leading to the conclusion that this system was very 
promising to study microlensing. Indeed, the predicted time de- 



* Based on observations made with the ESO-VLT Unit Telescope 
# 2 Kueyen (Cerro Paranal, Chile; Proposals 073.B-0243(A&B), 
074.B-0270(A), 075.B-0350(A), 076.B-0197(A), 177.B-0615(A&B), 
PI: F. Courbin). 



lays between the four quasar images are of the order of a day 
(Rix et al. 1992, Wambsganss & Paczynski 1994), meaning that 
intrinsic variability of the quasar can easily be distinguished 
from microlensing events. In addition, the particularly small 
redshift of the lensing galaxy implies large tangential veloci- 
ties for the microlenses. Furthermore the quasar images form 
right in the bulge of the lens where the stellar density is the 
highest. The combination of these properties makes microlens- 
ing events very likely in the Einstein Cross and very rapid, with 
timescales of a few weeks to a few months. Indeed, Irwin et al. 
(1989) reported significant brightness variations of the bright- 
est quasar image A, which they interpreted as the first detection 
ever of microlensing in the images of a multiply-imaged quasar. 

Since then, microlensing events have been observed in sev- 
eral other gravitationally lensed quasars, and are expected to 
occur in virtually any quadruply lensed quasar (Witt et al. 
119951 ). Probably the most compelling examples of microlens- 
ing light curves are given by the Optical Gravitational Lensing 
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Experiment (OGLE) (Wozniak et al. |2000al Udalski et al. 
2006). Started in 1997, this project monitors regularly the 
four quasar images of QSO 2237+0305, showing continuous 
microlensing-induced variations in the light curves. 

Most of the quasar microlensing studies so far are based ex- 
clusively on broad-band photometric monitoring (e.g. Wozniak 
et al. I2000bl Schechter et al. [20031 Colley & Schild 120031 . 
These observations, even though dominated by variations of 
the continuum, make it very difficult to disentangle variations 
in the continuum from variations in the broad emission lines 
(BELs). Both types of regions are affected by microlensing, 
but in different ways depending on their size. 

Microlensing of an extended source can occur when its size 
is smaller than or comparable to the Einstein radius of a star, 
i.e. of the order of 10 17 cm or 10~' pc in the case of the Einstein 
Cross (Nemiroff [T988l Schneider & Wambsganss |T990l ). From 
reverberation mapping, the broad line region (BLR) was long 
estimated to be larger than this, of the order of 10 18 cm or 
1 pc, hence leaving little room for BEL microlensing. However, 
more recent reverberation mapping studies revise this down- 
wards, to 10 16 cm (Wandel et al. [19991 Kaspi et al. 120001 
which is also consistent with the disk-wind model of Murray 
et al. d!9951 l. Inspired by these numbers, Abajas et al. (2002) 
and Lewis & Ibata (120041 ) investigated BEL microlensing in 
further detail and computed possible line profile variations for 
various BLR models. 

Observations of significant continuum and BEL mi- 
crolensing have been reported in a number of systems 
(QSO 2237+0305, Filippenko [19891 Lewis et al. [19981 
Wayth et al. [20031 HE 2149-2745, Burud et al. I2002al 
HE 0435-1223, Wisotzki et al. 120031 H 1413+117, Chartas 
et al. 120041 SDSS J1004+4112, Richards et al. 120041 
HE 1104-1805, Gomez-Alvarez 120041 HE 0047-1756, 
Wisotz ki et al. [20041 SDS S J0924+0219, Eigenbrod et al. 
I2006al Keeton et al. [20061 ; and RXJ 1131-1231, Sluse et 
al. 120071 ). These first observational indications of microlensing 
can be turned into a powerful tool to probe the inner parts of 
quasars, provided regular spectroscopic data can be obtained. 
Several theoretical studies show how multiwavelength light 
curves can constrain the energy profile of the quasar accre- 
tion disk and also the absolute sizes of the line-emitting re- 
gions (e.g., Agol & Krolik l 19991 Mineshige & Yonehara [l999l 
Abajas et al. l2002l Kochanek[2004]). 

In this paper we present the results of the first long-term 
spectrophotometric monitoring of the Einstein Cross. The spec- 
tral variations of the four quasar images are followed under 
sub-arcsecond seeing conditions with the Very Large Telescope 
(VLT) of the European Southern Observatory for more than 
two years, from October 2004 to December 2006, with a mean 
temporal sampling of about one point every second week. This 
first paper describes the observations, the method used to sepa- 
rate the quasar spectra from that of the lensing galaxy, and the 
main observational results. Simple considerations of the prop- 
erties of microlensing caustics and of the geometry of the cen- 
tral parts of the quasar already allow us to infer interesting con- 
straints on the quasar energy profile in the Einstein Cross. 

The full analysis of our monitoring data, still being ac- 
quired at the VLT, requires detailed microlensing simulations 



coupled with quasar models and will be the subject of future 
papers. 
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Fig. 1. VLT/FORS1 field of view showing the lensed quasar 
QSO 2237+0305, along with the four PSF stars used to spa- 
tially deconvolve the spectra. These stars are also used to cross- 
calibrate the observations in flux from one epoch to another and 
to minimize the effect of sky transparency. 



2. Observations 

We acquired our observations with the FOcal Reducer and low 
dispersion Spectrograph (FORS1), mounted on Kueyen, the 
Unit Telescope # 2 of the ESO Very Large Telescope (VLT) 
located at Cerro Paranal (Chile). We performed our observa- 
tions in the multi-object spectroscopy (MOS) mode. This strat- 
egy allowed us to get simultaneous observations of the main 
target and of four stars used as reference point-spread func- 
tions (PSFs). These stars were used to spatially deconvolve the 
spectra, as well as to perform accurate flux calibration of the 
target spectra from one epoch to another. We chose these stars 
to be located as close as possible to QSO 2237+0305 and to 
have similar apparent magnitudes as the quasar images. The 
PSF stars # 1, 2, 3, and 4 have R-band magnitudes of 17.5, 
17.0, 15.5, and 17.5 mag, respectively. Fig.[T]shows the field of 
view of our observations. 

We used the high-resolution collimator of FORS1 to 
achieve the best possible spatial sampling of the data, i.e. 0.1" 
per pixel. With this resolution, we observed a maximum of 8 
objects simultaneously over a field of view of 3.4' x 3.4'. One 
slit was aligned along two of the quasar images and four slits 
were centered on foreground PSF stars. We placed the remain- 
ing slits on empty sky regions and used them to carry out sky 
subtraction of the quasar data. 
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Fig. 2. FORS1 R-band acquisition image of QSO 2237+0305 
taken on epoch # 12 (12-09-2005). The slits used in the two 
masks are shown. Their width is 0.7". The Position Angle (PA) 
of mask 1 is PA = +56.5° and that of mask 2 is PA = +68.5°. 



Two masks were designed to observe the two pairs of 
quasar images. The PSF stars in both masks were the same. 
Fig.|2]shows the slit positioning with respect to our target. The 
first mask was aligned on quasar images A and D, while the 
second was aligned on images B and C. The masks were ro- 
tated to position angles that avoid clipping of any quasar im- 
age. This is mandatory to carry out spatial deconvolution of the 
spectra. 

Our observing sequences consisted of a short acquisition 
image, an "image-through-slit" check, followed by a consec- 
utive deep spectroscopic exposure. All individual exposures 
were 1620 s long. We list the journal of our observations in 
Table Q] The mean seeing during the three observing seasons 
was 0.8". We chose a slit width of 0.7", approximately match- 
ing the seeing and much smaller than the mean separation of 
1.4"between the quasar images. This is mandatory to avoid 
contamination of an image by the others. 

We used the G300V grism in combination with the GG375 
order sorting filter. For our slit width, the spectral resolution 
was AA — 15 A, as measured from the FWHM of night-sky 
emission lines, and the resolving power was R = A/AA 400 
at the central wavelength A = 5900 A. The useful wavelength 
range was 3900 < A < 8200 A with a scale of 2.69 A per pixel 
in the spectral direction. This configuration favors spectral cov- 
erage rather than spectral resolution, allowing us to follow the 
continuum over a broad spectral range, starting with the very 
blue portion of the optical spectrum. Even so and in spite of 
R = 400, a detailed profile of the BEL is still accessible. 

We also observed spectrophotometric standard stars 
(GD 108, HD 49798, LTT 377, LTT 1020, LTT 1788, and 
LTT 7987) to remove the response of the telescope, CCD, and 
grism. We did the relative calibration between the epochs using 
the PSF stars (see for more details Eigenbrod et al. 2006a). 



Table 1. Journal of the observations taken on 31 epochs. 



ID Civil Date HJD Mask Seeing ["] Airmass 



1 


13-10-2004 


3292 


1 


0.86 


1.204 


1 


14-10-2004 


3293 


2 


0.87 


1.221 


2 


14-11-2004 


3324 


1 


0.75 


1.184 


2 


14-11-2004 


3324 


2 


0.68 


1.305 


3 


01-12-2004 


3341 


1 


0.88 


1.355 


3 


01-12-2004 


3341 


2 


0.94 


1.609 


4 


15-12-2004 


3355 


1 


0.99 


1.712 


4 


16-12-2004 


3356 


2 


0.90 


1.817 


5 


11-05-2005 


3502 


1 


0.87 


1.568 


5 


12-05-2005 


3503 


2 


0.51 


1.389 


6 


01-06-2005 


3523 


1 


0.63 


1.342 


6 


01-06-2005 


3523 


2 


0.64 


1.224 


7 


01-07-2005 


3553 


1 


0.57 


1.153 


8 


14-07-2005 


3566 


1 


0.89 


1.620 


9 


06-08-2005 


3589 


1 


0.51 


1.135 


9 


06-08-2005 


3589 


2 


0.61 


1.173 


10 


15-08-2005 


3598 


1 


0.86 


1.140 


11 


25-08-2005 


3608 


1 


0.49 


1.261 


11 


25-08-2005 


3608 


2 


0.54 


1.461 


12 


12-09-2005 


3626 


1 


0.70 


1.535 


12 


12-09-2005 


3626 


2 


0.69 


1.341 


13 


27-09-2005 


3641 


1 


0.92 


1.480 


13 


27-09-2005 


3641 


2 


0.73 


1.281 


14 


01-10-2005 


3645 


1 


0.78 


1.281 


14 


01-10-2005 


3645 


2 


0.87 


1.156 


15 


11-10-2005 


3655 


1 


0.57 


1.140 


15 


11-10-2005 


3655 


2 


0.66 


1.134 


16 


21-10-2005 


3665 


1 


0.70 


1.215 


16 


21-10-2005 


3665 


2 


0.74 


1.156 


17 


11-11-2005 


3686 


1 


0.90 


1.137 


17 


11-11-2005 


3686 


2 


0.90 


1.185 


18 


24-11-2005 


3699 


1 


0.78 


1.265 


18 


24-11-2005 


3699 


2 


0.90 


1.443 


19 


06-12-2005 


3711 


1 


1.10 


1.720 


19 


06-12-2005 


3711 


2 


1.09 


1.445 


20 


24-05-2006 


3880 


1 


0.87 


1.709 


20 


24-05-2006 


3880 


2 


0.90 


1.443 


21 


16-06-2006 


3903 


1 


0.66 


1.213 


21 


16-06-2006 


3903 


2 


0.51 


1.155 


22 


20-06-2006 


3907 


1 


0.64 


1.286 


22 


20-06-2006 


3907 


2 


0.58 


1.189 


23 


27-06-2006 


3914 


1 


0.41 


1.145 


23 


27-06-2006 


3914 


2 


0.50 


1.133 


24 


27-07-2006 


3944 


1 


0.74 


1.316 


24 


27-07-2006 


3944 


2 


0.76 


1.204 


25 


03-08-2006 


3951 


1 


0.73 


1.246 


25 


03-08-2006 


3951 


2 


0.65 


1.169 


26 


13-10-2006 


4022 


1 


0.59 


1.176 


26 


13-10-2006 


4022 


2 


0.52 


1.300 


27 


28-10-2006 


4037 


1 


0.57 


1.148 


27 


28-10-2006 


4037 


2 


0.53 


1.138 


28 


10-11-2006 


4050 


1 


0.89 


1.515 


28 


10-11-2006 


4050 


2 


0.88 


1.323 


29 


27-11-2006 


4067 


1 


0.87 


1.255 


29 


27-11-2006 


4067 


2 


0.92 


1.391 


30 


19-12-2006 


4089 


2 


1.04 


2.125 


31 


22-12-2006 


4092 


1 


0.80 


2.018 


31 


23-12-2006 


4093 


2 


0.76 


2.248 



4 



A. Eigenbrod et al.: Microlensing variability in the Einstein Cross 




Fig. 3. Spectra of the four PSF stars. The spectra in each panel 
correspond to different observing epochs, chosen to span the 
full length of the monitoring. The IDs of the observing epochs, 
as given in Table [T] are indicated. The differences in flux are 
mainly due to the presence of thin clouds. The purpose of us- 
ing these stars as flux cross-calibrators is precisely to eliminate 
these differences, both in intensity and shape. 



Finally, it is worth emphasizing that all our VLT data used 
in the present paper were obtained in service mode, without 
which this project would have been impossible. 



Fig. 4. Flux correction for different epochs with respect to the 
reference epoch # 2 (14-11-2004). Each panel corresponds 
to one of the 4 PSF stars visible in Fig. [3] In each panel, the 
dotted line shows the ratio of the spectrum of one PSF star 
taken at a given epoch and the spectrum of the same star taken 
at the reference epoch # 2. The curves are polynomials fitting 
the data. Importantly, the correction derived at a given epoch is 
about the same for the four stars. The mean of these four curves 
is used to correct the spectra of the Einstein Cross images. The 
small parts of the spectra with strong atmospheric absorption 
are masked. The different spectral ranges are due to different 
clippings of the spectra by the edges of the CCD. 



3. Data analysis 

3.1. Reduction 

The data reduction followed the same procedure described in 
detail in Eigenbrod et al. (2006b). We carried out the standard 
bias subtraction and flat field correction of the spectra using 
IRAF3. We obtained the wavelength calibration from the spec- 
trum of helium-argon lamps. All spectra, for the object and for 
the PSF stars were calibrated in two dimensions. 

Only one single exposure was taken per mask and per 
epoch. For this reason, the usual cosmic-ray rejection scheme 
applied to multiple images could not be applied. Instead, we 
used the L. A. Cosmic algorithm (van Dokkum |2001| >, that 
can handle single images. We visually inspected the cosmic- 
ray corrected images to check that no data pixel was affected 
by the process, especially in the emission lines and in the data 
with the best seeing. 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



We removed the sky background in a different way in the 
spectra of the PSF stars and in those of the gravitational lens. 
For the PSF stars, which are small compared with the slit length 
(19"), we used the IRAF task background. This task fits a 
second order Chebyshev polynomial in the spatial direction to 
the areas of the spectrum that are not illuminated by the ob- 
ject, and subtracts it from the data. As the lensing galaxy in 
QSO 2237+0305 is larger than the slit length, this procedure is 
not applicable. Instead, we used the slits positioned on empty 
sky regions of the FORS 1 field of view, and located next to the 
gravitational lens. The sky was fitted to these slits and removed 
from the slit containing the images of QSO 2237+0305. 

3.2. Flux cross-calibration 

Once the cosmic rays and sky background were removed, we 
applied a flux cross-calibration of the spectra as described by 
Eigenbrod et al. (2006a), using the four PSF stars. The spectra 
of these stars are shown in Fig. [3] for five different observing 
epochs. Our observations show that these stars are non variable. 



A. Eigenbrod et al.: Microlensing variability in the Einstein Cross 



5 




Fig. 7. Deconvolved and extracted ID-spectra of the lensing 
galaxy. The two panels correspond to the two MOS masks. 
The shaded areas are the envelopes containing all the spectra 
of the lens obtained with the corresponding mask. The thick 
black lines are the means. Note the small scatter between the 
two spectra. 



Fig. 8. Example of a spectral decomposition. The top panel 
shows the two extracted spectra for the images A and D for 
the observations taken on epoch # 25 (03-08-2006) with mask 
1 . The extracted spectrum of the lensing galaxy, in the bottom 
panel, shows no trace of contamination by the quasar BELs. 
For clarity, we have multiplied the flux of image D by a factor 
of two. 



We created a ratio spectrum for each star, i.e. we divided the 
spectrum of the star by the spectrum of the same star for a cho- 
sen reference exposure. We chose epoch # 2 (14-1 1-2004) as 
our reference exposure because of the excellent weather con- 
ditions at this particular epoch for both seeing and sky trans- 
parency. The computation of these flux ratios was done for all 
four stars in each exposure and we checked the compatibility 
of the response curves derived with the four different stars are 
compatible (see Fig.|4]i. If not, we rejected one or a maximum 
of two of the PSF stars. This can happen in some exceptional 
cases, e.g. when the alignment between the star and the slit is 
not optimal and generates a color gradient in the spectrum of 
the misaligned object. Aside from this instrumental effect, the 
observations show no trace of intrinsic variability of the PSF 
stars. 

After checking that the correction spectra obtained for the 
four stars were very similar, we computed their mean, which 
we took as the correction to be applied to the gravitational lens. 
The high stability of the corrections across the field demon- 
strates that all residual chromatic slit losses due to the atmo- 
spheric refraction are fully corrected. This correction is eased 
by the fact that: (1) the position angle of the masks is the same 
for the quasar images and for the PSF stars (i.e. the PSF clip- 
ping is the same for the target and the reference stars); (2) 
we avoid observations at large airmasses (i.e. never larger than 
2.5); and (3) the atmospheric refraction corrector on FORS1 is 
very efficient. 



3.3. Deconvolution 

The lensing galaxy in QSO 2237+0305 is bright. Its 
central parts have a surface brightness of approximately 
18 mag/arcsec 2 in the R band, which is comparable to the 
quasar images. Hence studying microlensing variations of the 
quasar images requires very accurate deblending. 

In order to carry out this challenging task, we used the spec- 
tral version of the MCS deconvolution algorithm (Magain et al. 
1998, Courbin et al. 2000), which uses the spatial information 
contained in the spectra of the PSF stars. The algorithm sharp- 
ens the spectra in the spatial direction, and also decomposes 
them into a "point-source channel" containing the spectra of the 
two quasar images, and an "extended channel" containing the 
spectrum of everything in the image that is not a point source, 
in this case, the spectrum of the lensing galaxy. In Fig. [5] we il- 
lustrate the process and the different outputs. In Fig. [7] we show 
how similar the spectra of the lensing galaxy are, extracted ei- 
ther from two different masks or from data taken at different 
epochs, hence illustrating the robustness of the deconvolution 
technique. In Fig. [8] we give an example of decomposition of 
the data into the quasar and lens spectra after integrating along 
the spatial direction. The lensing galaxy spectrum shows no 
trace of the residual quasar BELs. Even when the contrast be- 
tween the quasar and the galaxy is particularly large, the de- 
composition is accurate. For example, the Call H+K doublet in 
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Fig. 5. Left: portion of the VLT 2D-spectrum of quasar images D and A, taken on epoch # 25 (03-08-2006), on which are 
indicated the main spectral features of either the quasars or the lens. Center left: spatially deconvolved spectrum. The two quasar 
images are very well separated. Center right: spectrum of the lensing galaxy alone. Right: residual map of the deconvolution 
after subtraction of the quasar and lens spectra. Note that the residuals are displayed with much narrower cuts than those used 
in the other panels. The darkest and brightest pixels correspond to -3cr and +3cr respectively. No significant residuals of the 
spectral features are visible. 




U Lj i i I : : i i i i I i i i i I i i LJ i I I 

4000 5000 6000 7000 8000 

A [A] 

Fig. 6. Deconvolved and extracted spectra of quasar image A for five observing epochs. Chromatic variations in the spectra are 
conspicuous with the blue part of the spectra being more magnified than the red part. 



the lens spectrum is well visible, in spite of the presence of the 3.4. Cross-check with the OGLE-HI light curves 

strong quasar C IV emission in the same wavelength range. 

After reduction and spatial deconvolution, we obtained the 
extracted spectra of quasar images A and D on 30 different 
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Fig. 9. OGLE-III light curves (Udalski et al. [2006b of all four 
quasar images from April 2004 to December 2006 (dots), 
compared with the photometry derived by integrating our 
VLT spectra through the OGLE V-band (dark triangles). The 
1-sigma error bars correspond to the photon noise in the 
spectrum. We shift the OGLE-III light curve of image D by 
-0.5 mag with respect to the published values. The bottom 
panel displays the seeing values for each observations. 



epochs, and of B and C on 28 different epochs. Several ex- 
tracted spectra of image A are shown in Fig|6] As a sanity 
check, we compared our results with the OGLE-III photomet- 
ric monitoring of QSO 2237+0305 (Udalski et al. 120061 We 
integrated our quasar spectra in the corresponding V-band to 
estimate, from the spectra, the photometric light curves as if 
they were obtained from imaging. In Fig. [9] we compare our 
magnitude estimates with the actual OGLE-III measurements. 
The overall agreement is very good for images A, B, and C. 
For image D, we have to shift the OGLE-III light curve by 
-0.5 mag with respect to the published values. Interestingly, 
this shift is not needed when we compare our results with the 
previous OGLE data from the provisional calibration presented 
in the years 2004-2006. The previous OGLE data also agreed 
with the photometry of Koptelova et al. (|2005). This changed 
when Udalski et al. (2006) reviewed their calibration and gave 
image D a larger magnitude of approximately 0.5 mag. They 
stated that the steep rise of brightness of image D at the end 
of the 2000 OGLE-II season leaded to an overestimate of the 
extrapolated magnitude for the beginning of the 2001 OGLE- 
III season. But this is now discrepant with the photometry of 
Koptelova et al. (120051 1 . We think that the new extrapolation of 
the light curve of image D from the end of season 2000 to the 
beginning of season 2001 might be uncertain, leading to the 
observed shift between our data and the OGLE-III light curve 



of image D. However, aside from this shift, the agreement be- 
tween the OGLE photometry and our integrated VLT spectra is 
also very good for image D. 

4. Multi-component decomposition 

Different emission features are known to be produced in re- 
gions of different characteristic sizes. As microlensing magni- 
fication varies on short spatial scales, sources of different sizes 
are magnified by differing amounts (e.g. Wambsganss et al. 
119901 1. Emission features from smaller regions of the source are 
more highly variable due to microlensing than features emitted 
in more extended regions. In order to study the variation of 
each spectral feature independently, we need to decompose the 
spectra into their individual components. 

4.1. Method 

In our analysis of the 1-D spectra of the four quasar images, we 
follow the multi-component decomposition (MCD) approach 
(Wills et al. [19851 Dietrich et al. 120531 implemented in Sluse 
et al. (2007). This method is applied to the rest-frame spectra, 
assuming they are the superposition of (1) a power law con- 
tinuum, (2) a pseudo-continuum due to the merging of Fe II 
and Fe III emission blends, and (3) an emission spectrum due 
to the other individual BELs. We consider the following emis- 
sion lines : C IV A 1 549, He II Al 640, III] A 1 664, Al II A 1 67 1 , 
A1III /H857, Si III] 1892, CIH] ,11909, and Mgll ,12798. All 
these features are fitted simultaneously to the data using a stan- 
dard least-square minimization with a Levenberg-Marquardt 
based algorithm adapted from the Numerical Recipes (Press 
etal. [19861 

In the first step, we identify the underlying nonstellar 
power-law continuum from spectral windows that are free 
(or almost free) of contributions from the other components, 
namely the iron pseudo-continuum and the BELs. We use the 
windows 1680 < A < 1710 A and 3020 < A < 3080 A. 
After visual inspection of the iron templates by Vestergaard et 
al. ([2001 >, we do not expect significant iron emission in these 
windows. 

We characterize the spectral continuum (measured in the 
restframe) with a power law /„ oc V* v , which translates in wave- 
length to fx A ax with the relation a v = -(2 + o^), i.e. 

where Aq = 2000 A and where a v is the commonly used canon- 
ical power index. 

Next, we fit the BELs with Gaussian profiles. We consider 
a sum of three profiles to fit the absorption feature in the C IV 
emission line. Two profiles are used for the C III] line and one 
single profile is used to fit simultaneously the O III] and Al II 
lines. All other BELs are fitted with one single profile. We 
then subtract the BELs and the continuum from the spectra. 
We consider the residuals as coming from the emission blends 
of Fe II and Fe III. Hence the averaged and normalized residu- 
als over all epochs define our first iron pseudo-continuum tem- 
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Fig. 10. Multi-component decomposition of the spectrum of the brightest quasar image, A, taken on epoch # 2 (14-1 1-2004). 
The upper panels show the detailed spectral decomposition of the BELs, while the middle panel displays the entire spectrum. 
The continuum is indicated as a dotted curve. The Gaussian lines and iron pseudo-continuum templates are shown below the 
spectrum. The bottom panel is the residual for each pixel normalized by the photon noise per pixel (i.e. the y-axis is the residual 
flux in units of <x). 




Fig. 11. Left: Examples of light curves for the quasar image A. The integrated flux for the continuum, the He II, and C III] BELs 
are given from top to bottom. The continuum is integrated over the entire available wavelength range. In each panel, we fit a scaled 
version (solid line) of the OGLE-III light curve (Udalski et al. 2006). This nicely illustrates that the BELs vary simultaneously 
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Fig. 12. Correlation between the intensity fo and slope a v of 
the continuum spectra for all four quasar images. The points 
are connected chronologically. The first observation epoch is 
marked by a star and the last one by a square. The correla- 
tion is obvious in images A and B spanning a broad range of 
spectral slopes. In these two images, an increase in intensity 
is accompanied by an increase in steepness, i.e. when a quasar 
image gets brighter, it also gets bluer. 



images A and B, indicating that significant microlensing events 
occured in these two images. 

Inaccurate alignment of the quasar images in the slit of the 
spectrograph is a possible instrumental effect that can mimic 
microlensing changes in the spectral slope of the quasar im- 
ages. Indeed, small clipping of one of the quasar images would 
lead to a stronger flux loss at the bluer wavelengths, hence 
producing a color gradient in the spectrum and a decrease in 
the measured value a v . We have checked all the "through-slit" 
images taken before each spectrum. Not only do these images 
show that the alignment is correct, but it is also very easily re- 
producible from one epoch to another, even when the FORS 1 
has been dismounted from and remounted on the telescope. 

We have also checked that our fitting procedure does not in- 
troduce any spurious correlation between a v and fo- We check 
this by using simulated spectra. In order to do that we take a 
reference spectrum for each quasar image and subtract its con- 
tinuum. We then take random pairs of (a v , fo) parameters so 
that the a v vs. fo plane is well sampled. We chose 400 such 
pairs and create the corresponding continuum to be added to 
the reference spectrum. The decomposition procedure is then 
run on the 400 spectra. We find no correlation at all between 
the measured a v and fo. In addition, the parameters used to 
build the simulated spectra are almost perfectly recovered by 
the decomposition procedure. 

We conclude that genuine chromatic variations are present 
in the continuum of all images of QSO 2237+0305. The effect 
is most pronounced in image A during the last observing sea- 
son, and in image B at the beginning of our monitoring. We 
show in the following that these observed variations are, in ad- 
dition, well compatible with the OGLE-III single-band photo- 
metric observations. 



plate. We can then proceed iteratively by including this pseudo- 
continuum iron template in the fitting procedure and rerun it. 
This gives a better fitting of the emission lines and defines an 
improved iron pseudo-continuum template. After five of these 
iterations, typically the fitting does not change significantly 
anymore. Fig. [10] shows an example of the fitting decompo- 
sition described above. 

4.2. Results 

The light curves for the continuum and for the emission lines 
can be constructed from the above multi-component decom- 
position. We show in Fig. [TT] an example of variation in the 
brightest quasar image, A, for the continuum and for two BELs. 
The error bars give the photon noise, integrated over the corre- 
sponding wavelength range. In the right panel of Fig. [TT] we 
show the variability of the continuum in intensity, fo, and in 
slope, a v , for the 4 quasar images. It is immediately clear that 
the continuum variations with the largest amplitude are ob- 
served in image A, between HJD = 3600 and 3900 days, and 
in image B between HJD = 3300 and 3500 days. These varia- 
tions are accompanied by an increase in steepness, i.e. when a 
quasar image gets brighter, it also gets bluer. This is particularly 
obvious in Fig.Q~2] where fo and a v are strongly correlated for 



5. Microlensing variability in the OGLE-III 
photometry 

The photometric variations in most gravitationally-lensed 
quasars are dominated by the intrinsic variations of the quasar, 
typically of the order of 0.5 - 1.5 mag, hence making them 
useful to measure the time delays between the quasar images. 
Microlensing variations are usually smaller, in the range 0.05 - 
0.1 mag (e.g. the lensed quasars B 1600+434 by Burud et al. 
120001 RX J091 1+0551 by Hjorth et al. l2002l SBS 1520+530 
by Burud et al. I2002bl FBQ 0951+2635 by Jakobsson et al. 
120051 

The Einstein Cross is different from this general behavior 
in two ways: (1) the time delays between each pair of images 
are expected to be of the order of one day, hardly measurable, 
and (2) the microlensing variations dominate the light curves. 
For these two reasons, microlensing can be fairly well isolated 
in each quasar image, because it acts differently on the four 
sightlines. 

To separate the intrinsic flux variations of the quasar from 
the microlensing ones, we perform a polynomial fit to the 
OGLE-III light curves (Udalski et al. 120061 of Fig. [9] This 
simple and fully analytical method has been developed by 
Kochanek et al. (2006|, and is also described by Vuissoz et al. 
(2007). In the present application, the variations of each quasar 
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Fig. 13. Decomposition of the OGLE-III photometric light 
curves (Udalski et al. [2006) of the quasar images, into intrin- 
sic quasar variations and microlensing-induced variations (see 
Section|5}. The intrinsic variations are shown at the bottom of 
the figure as a continuous line, while the pure microlensing 
variations are the data points. The curves are shifted arbitrar- 
ily along the y-axis for clarity. The tickmarks at the top show 
the epochs of our observations. 



image are modeled as a sum of two Legendre polynomials: one 
polynomial is common to all four quasar images and represents 
the intrinsic variations of the source while a second polyno- 
mial, different for each quasar image, represents the additional 
microlensing variations. In doing so, we rescale the OGLE-III 
error bars of each image by a factor equal to the flux ratio be- 
tween each image and image C. This rescaling suppresses the 
potential problem existing if the fitting procedure considers the 
variation of image A (with the highest signal-to-noise) as the 
intrinsic variation of the quasar. The chosen order of the poly- 
nomial is 7 for the intrinsic variation, and 10 for the microlens- 
ing variation. Higher order polynomials do not significantly im- 
prove the fit. The results are displayed in Fig. [13] where the in- 
trinsic variation of the source recovered by the simultaneous fit 
is shown together with the pure microlensing variations. 

We check the efficiency of our method by generating arti- 
ficial light curves and then using the above polynomial fit to 
recover the intrinsic and microlensing light curves. These arti- 
ficial light curves are generated in the same way as described 
in Eigenbrod et al. (120051 1. and are composed of an intrinsic 
light curve to which we add microlensing fluctuations. Both 
are created in a random walk manner (i.e. not from polynomi- 
als). They are constructed to match the variability properties 
of the actual light curves, i.e. their timescale and amplitude of 
variation (for further details see Eigenbrod et al. 2005). We re- 



cover the simulated intrinsic light curve of the quasar with a 
typical error of less than 0.1 mag. The variations of more than 

0. 4 mag, shown in Fig. [13] both for microlensing and quasar 
variations, are well above the error estimated from the simu- 
lated light curves. In our simulations, we adopt the same pho- 
tometric error bars as in the light curves of all quasar images, 

1. e. the re-scaling of the error bars described above in the real 
data is taken into account. If, on the contrary, we adopt error 
bars that follow the photon noise, the fitting procedure consid- 
ers the highest signal-to-noise light curve as the intrinsic quasar 
light curve. 

The light curve most affected by microlensing is that of 
image B, with a peak-to-peak amplitude of more than 0.7 mag 
over 3 years. The other quasar images show microlensing- 
induced variations of up to 0.4 mag, with quasar image A hav- 
ing a sharp event during the last observing season. The intrinsic 
quasar light curve displays a variation of about 0.4 mag. 

The polynomial decomposition of the light curves are com- 
patible with the spectroscopic results. Quasar images A and B, 
which have the largest microlensing contribution in Fig. [13] re- 
spectively at HJD ~ 3500 days and HJD ~ 3900 days, also have 
a sharp rise in a v at the same epochs. 

6. Microlensing variability in the spectra 

Chromatic variations of the continuum of images A and B of 
QSO 2237+0305 are clearly seen in our data. In addition, dif- 
ferential magnification of the continuum with respect to the 
BELs is also seen in all four quasar images. Such effects have 
already been observed by Lewis et al. ( 1998) and Wayth et al. 
( 2005 ), but only for data over two epochs. Our VLT spectra al- 
low us to follow the variations over two full years, provided the 
intrinsic variations of the quasar are removed. 

6.1. Continuum and BELs relative magnifications 

Since the time delays in QSO 2237+0305 are negligible, taking 
the ratio between the above quantities in pairs of quasar images 
cancels the intrinsic variations. Let F(t) be the intrinsic source 
flux, and M,, yu, be the macro and microlensing-magnifications 
of quasar image i, respectively. The observed flux ratio between 
images i and j at time t is then : 



Rij(t) 



Hi(f) Mj Fit) 
Hj(t) Mj F{t) 



W(t) Mj 
fij(t) Mj 



<r t -ri) 



(1) 



The extinction e~ Ti remains constant in time and is rela- 
tively similar in all four quasar images as we show in Sec. 16.31 
Hence, it is not expected to strongly affect our results and we 
will neglect it in the following. The macromagnifications M, 
are best estimated in the mid-IR and radio domain (Falco et al. 
1996, Agol et al. 2000). At these wavelengths, the source size 
is much larger than the typical spatial scale in the microcaus- 
tics network, hence leaving it fairly unaffected by microlensing 
(i.e. fii = 1). By multiplying Rjj by Mj/Mi, using the mid-IR 
observations, we find the pure-microlensing magnification ra- 
tios: 



r tj (t) = 



(2) 
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Fig. 14. Microlensing-magnification ratios ry as a function of time (HJD-2450000) for all possible combinations of image pairs. 
These ratios are corrected for the macromagnification by the lensing galaxy following equation ©. The ratios are given for the 
continuum, for the BELs He II, C IV, C III], Mg II and for the iron pseudo-continuum. 
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In Fig. [14] we show the variations of r, ; (f) for the integrated 
flux in the main emission lines and in the continuum. We plot 
in the figure all possible combinations of image ratios. Several 
interesting results can be drawn from this figure: 

(a) With the exception of r BD during the last season (i.e. 
HJD > 3900 days), none of the ratios r, ; are close to 
1, meaning that microlensing is acting at least in three 
quasar images during the entire monitoring. The evidence 
for rgD ~ 1 when HJD > 3900 days supports the absence of 
significant microlensing in images B and D during the last 
observing season. 

(b) All the BELs have + 1 (except r^oi- The r, ; ratios in 
the BELs generally closely follow the value measured for 
the continuum. This demonstrates that the BLR is small 
enough - probably not larger than a few Einstein radii of 
a typical microlens in QSO 2237+0305 - so that the BELs 
can be significantly affected. In addition, the variations ob- 
served in the BELs are correlated with those in the contin- 
uum. 

(c) The largest changes of magnification ratios involve images 
A (for HJD ~ 3900 days) and B (for HJD ~ 3500 days). 
This is seen already in Fig. [13] but the ratios shown here 
allow us to make sure that the variations are only due to 
microlensing. 

(d) The ratios r&j and tq deviate significantly from 1 during 
the whole observing period. As an example, we measure 
in the continuum (resp. CHI]) r&D ~ 2.6 (resp. 2.8) for 
HJD < 3900 days and rcD ~ 1.6 (resp. 1.6). Interestingly, 
these ratios have remained roughly the same since the be- 
ginning of the OGLE-HI campaign (i.e. HJD ~ 2000 days). 
Our values of tad and tqd are also similar to the value mea- 
sured in the H /3 emission line in August 2002 (HJD = 2496 
days) by Metcalf et al. (2004). This strongly suggests that 
images A and C are affected by long-term micro/milli- 
lensing on periods longer than 5 years. 

(e) The ratio rcD is the most stable ratio along the monitor- 
ing campaign, indicating that no major (short) microlensing 
event occured in images C or D. 



Table 2. Mean microlensing ratios for the continuum and for 
the main BELs. The mean values for rAD = haIhd and rgD = 
UbIhd are computed for the observations around the epoch in 
the HJD line, i.e. during microlensing events or during quies- 
cent phases. The values are given along with the dispersion of 
the points around the mean. 





< I'AD > 


< r AD > 


< t'BD > 


< f ED > 


HJD 


3900 d 


3500 d 


3500 d 


3300 d 


State 


Micro-A 


Quiet-A 


Micro-B 


Quiet-B 


Cont. 


3.46 ± 0.24 


2.65 ±0.18 


1.28 ± 0.12 


0.79 ± 0.05 


CIV 


2.89 ±0.13 


2.90 ±0.11 


1.22 ± 0.11 


1.15 ±0.07 


Hell 


3.01 ± 0.55 


2.70 ± 0.33 


1.50 ±0.23 


1.39 ±0.10 


CIII] 


2.71 ± 0.07 


2.54 ± 0.09 


1.17 ± 0.09 


1.02 ± 0.04 


Mgll 


2.49 ± 0.10 


2.40 ±0.10 


0.88 ± 0.09 


0.91 ± 0.04 


Fell+III 


2.56 ± 0.23 


2.05 ±0.17 


0.99 ±0.15 


0.75 ± 0.09 



ionization potential (i.e. Mg II). Indeed the ionization poten- 
tials of the different lines are 47.9 eV (CIV), 24.6 eV (Hell), 
24.4 eV (CIII]), and 7.6 eV (Mgll). The difference of magni- 
fication between C III] and Mg II was not observed by Wayth 
et al. d20051 >. For image B, the same global trend is observed 
except that the relative errors on the rgD ratios are higher due 
to the lower signal-to-noise ratio of the spectra of image B. The 
effect is almost absent when r^D and rgD are computed in qui- 
escent phases of components A and B. 

The behaviour of the Fell+III emission is more difficult 
to interpret, as this complex is in fact a blend of many lines. 
However, we note that the Fe II+III complex in image A is mi- 
crolensed at about the same level as the Mg II line. In addition, 
the difference in magnification of the Fe II+III lines between 
a microlensing and a quiet phase is larger than for the other 
lines. This may suggest differential magnification of the emit- 
ting regions within the Fe II+III complex, i.e. that the Fe II+III 
is present both in compact and more extended regions, a con- 
clusion also reached by Sluse et al. (12007b . 



Because of (b) and (e), we can safely consider that during 
the time span of our observations, image D is the less affected 
by microlensing, both on short (i.e. of the order of a few weeks) 
and long timescales (years). This is consistent with the broad- 
band microlensing light curves of Fig.fOl 

We use image D as a reference to study the "short" mi- 
crolensing events affecting image A at HJD ~ 3900 days and 
image B at HJD ~ 3500 days. During both events the emission 
lines are not magnified by the same amount f Fig. fL4b. 

To quantify this, we compute in Table [2] the mean values 
of rAD and rsD during microlensing events and in more "quies- 
cent" phases. For image A (i.e. tad), we clearly see in this table 
that C IV, C III] and He II show very similar magnification ra- 
tios, while the Mg II line is less magnified. 

All the lines are less magnified than the continuum, consis- 
tent with a scheme where the continuum is emitted in the most 
compact region, and where other emission lines are emitted in 
larger regions, the largest region being the one with the lowest 



6.2. Line profiles 

We have discussed in the previous sections the global inten- 
sity changes in the emission lines of QSO 2237+0305. We 
now investigate the possibility of a change in their profile. 
Such profile variations may be caused by differential magni- 
fication of regions with different velocities in the BLR. This 
can introduce, e.g., asymmetric changes or even peak dis- 
placement of the line (Lewis 119981 Abajas et al l20021 Lewis 
& Ibata 2004). Observational evidence for profile variations 
has been reported, e.g. in HE 2149-27 (Burud et al. I2002ab : 
SDSS 1004+4112 (Richards et al. 120041 Gomez-Alvarez et 
al. 120061 ; SDSS J0924+0219 (Eigenbrod et al. I2006at ; and 
RXJ 1131-1231 (Sluse et al. l2007l l. 

We study the variability of the emission lines in the four 
quasar images using the continuum-subtracted spectra obtained 
from the decomposition procedure of Sec. 14. II We concentrate 
on the C III] emission line, which is easily decomposed into 
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Fig. 15. The C IV (left) and C III] (right) BEL profiles for quasar images A, B, C and D as observed at epoch # 1 (14-10-2004, 
dotted) and epoch # 24 (27-07-2006, solid). The first epoch is used as a reference while the second epoch falls within the high- 
magnification episode that occured at HJD ~ 3900 days in image A (Sect. [6721 . In order to properly show the extra-magnification 
of the line wings during the strong microlensing episode of image A, we have multiplied the line profile in A observed on 
14-10-2004 by a factor 1.05. 



a sum of two components with full-width-half-maximum val- 
ues about 1 100 km s _I and 5 000 km s _1 . We will refer to these 
components in the following as the "narrow" and "broad" com- 
ponents, although the difference made here between narrow 
and broad is only phenomenological. In particular, the narrow 
component is not associated, a priori, with the narrow line re- 
gion (NLR). 

We measure the microlensing magnification ratio 
p(,(narrow)/yU,(broad) for the CHI] emission line in each 
quasar image. In image A, the magnification of the broad 
component is ~1.8 times larger than the magnification of the 
narrow component, and this does not change drastically along 
our whole monitoring campaign. In the other quasar images, 
the microlensing-magnification of the narrow and broad com- 
ponents are comparable. This implies that the C III] emission 
line is microlensed globally in all images, except in A. In the 
latter image, the broadest part of the line is more microlensed 
than the core, indicating that the broadest component of the 
C III] emission is emitted in a more compact region than the 
core of the emission line. 

The core of the narrow C III] emission line is probably as- 
sociated mainly with photons emitted by the NLR: we observe 
nearly the same magnification ratio between the narrow and the 
broad components of the C III] emission line than is reported 
between the narrow [OIII] and the broad H/3 emission lines 
(Metcalf et al. 2004). This gives a hint that the narrow compo- 
nent of the C III] emission line may be partly emitted by the 
NLR. 

However, we do not measure the same microlensing- 
magnification ratios ry in the narrow C III] emission as those 
measured by Metcalf et al. (2004) in the narrow [OIII] emis- 



sion line. The rgc, rsD, rcD ratios measured by these authors 
are similar for the H/3 BEL and for the narrow [OIII] emis- 
sion line, but are different from 1, which may be interpreted 
as a consequence of the microlensing of the NLR. If the NLR 
is indeed microlensed, the discrepancy between our values and 
the ones by Metcalf et al. (j2Q04S) can simply be explained by 
the fact that the amplitude of microlensing changes with time. 
However, there are also two other possible explanations. First, 
the size of the [OIII] emission region might be comparable 
with the size of the macrocaustic of the lensing galaxy, such 
that the [O III] lensed images are extended and can be resolved 
by integral field spectroscopy (Metcalf et al. 2004, Yonehara 
2006). This may lead to uncertainties in the flux measurements 
depending on the chosen size of the aperture. Finally, the dis- 
crepancy may be explained as well in terms of extinction by the 
lens (see Sec. 16.31 ). 

The details of the C III] profile also show variations. Our 
spectra show evidence for systematic broadening, by ~10% of 
the C III] emission line in quasar image A during our last ob- 
serving season (HJD > 3900 days; Fig. [T5l ). In addition, we 
note evidence for variations in the central parts of the CIV 
emission lines in all components and at most epochs. An ex- 
ample of these variations is shown in Fig. [15] Their interpreta- 
tion will be much more complex than for the C III] due to the 
decomposition of the line into three emission components plus 
an absorber-like feature. 

To look for distortions in the line profiles of images B, C, 
and D, we normalize the continuum-subtracted spectra of all 
images so that they share the same C III] central intensity. The 
result is shown for selected epochs in Fig. [16] With this choice 
of normalization and in the absence of microlensing, the C III] 
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Fig. 16. Comparison of the C III] BEL profiles of images A, B, C and D of the lensed quasar QSO 2237+0305 on 6 different 
epochs. The line profiles are normalized so that they have the same peak value. This representation illustrates the difference in 
the line profiles of the 4 lensed images at several epochs. The most striking effect is the larger wings observed at all epochs in 
image A. These larger wings are caused by differential microlensing between the wings and the core of the CIII] line. Small 
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lines at a single epoch would match perfectly. It is conspicuous 
in FigfT6]that this is not the case. First, we clearly see the effect 
of differential magnification between the core and the wings 
of image A at all epochs (the wings of C III] in image A are 
always larger than in the other emission lines). The emission 
lines in B, C, and D are more similar to one another but there is 
no epoch where the three line profiles match perfectly, indica- 
tive that small microlensing fluctuations affect the C III] line. 
Image C also shows line profile variations, even though the ef- 
fect is less pronounced than in image A. Finally, we note the 
absence of any strong line profile variations of image B dur- 
ing the short-term microlensing event occuring in that image at 
HJD ~ 3500 days. 

6.3. Differential extinction by the lensing galaxy 

Most of the magnification ratios fluctuate around a mean 
value. The fluctuations themselves can only be explained by 
microlensing, but the value of the mean < r, ; > is usually dif- 
ferent from 1 during our observations. This can, in principle, be 
explained either by long-term microlensing or by non variable 
extinction by dust in the lens. 

In order to test the latter hypothesis, we use a simple and 
empirical diagnostic using ratio spectra of pairs of quasar im- 
ages. In the absence of reddening and microlensing, these ratio 
spectra should be flat. If dust is present in different amounts on 
the lines of sight of the images, the ratio spectra will show a 
non-zero slope, constant with time. Any time-variable change 
of slope can safely be attributed to microlensing. 

We find that the most useful pairs are the ones formed by 
A&C, B&C, and C&D. Indeed, the C/D ratio spectrum is found 
to be almost flat all along the years, indicating no significant 
differential extinction. This is not the case for the two other 
pairs of ratio spectra which show reddening of image C relative 
to both A and B, as also reported by Yee d 1 988t . We estimate 
this differential reddening using the extinction law by Cardelli 
et al. (1989) and by assuming R v = 3.1. We find that the dif- 
ferential extinction A V (C) - A V {A) A V (C) - A V (B) is in the 
range 0.1-0.3 mag. This range of values is sufficient to explain 
the discrepancy found between the ry ratios measured in the 
C III] narrow component and in [O III] of Sec. 16.21 

Finally, our estimates of the extinction in the lens are too 
small to explain the highest values of the mean magnification 
ratios observed in Fig. [14] and Table [2] For instance, the mean 
of the tab or r^D ratio reaches values larger than 2. Static, long- 
term microlensing is therefore present in the Einstein Cross, at 
least in images A and C. 

7. Conclusions 

This paper presents the first long-term (2.2 years) and well- 
sampled spectrophotometric monitoring of a gravitationally 
lensed quasar, namely the Einstein Cross QSO 2237+0305. 
The mean temporal sampling is of one observation every sec- 
ond week. The observations are carried out with the VLT in a 
novel way, using the spectra of PSF stars, both to deblend the 
quasar images from the lensing galaxy and to carry out a very 



accurate flux calibration. This paper, the first of a series, de- 
scribes the observations and the techniques used to extract the 
scientific information from the data. 

Detailed inverse ray-shooting simulations will be needed to 
infer quantitative information on the internal structure of the 
lensed quasar, and will be the topic of the future papers. The 
main observational facts that these simulations will need to take 
into account can be summarized as follows. 

We find that all images of QSO 2237+0305 are af- 
fected by microlensing both on the long and short timescales. 
Comparison of the image flux ratios with mid-IR measure- 
ments reveals that quasar images A and C are both affected 
by long-term microlensing on a period longer than 5 years. 
This long-term microlensing affects both the continuum and 
the BELs. 

Furthermore, in quasar image A, the broad component of 
the C III] line is magnified by a factor 1.8 larger than the narrow 
component. On the contrary, the other quasar images have the 
same magnification in the narrow and broad components. 

On the short timescales, i.e. several months, images A and 
B are the most affected by microlensing during our monitor- 
ing campaign. Image C and especially D are the most qui- 
escent. Image A shows an important brightening episode at 
HJD ~ 3900 days, and image B at HJD ~ 3500 days. We show 
that the continuum of these two images becomes bluer as they 
get brighter, as expected from microlensing magnification of an 
accretion disk. 

We also report microlensing-induced variations of the 
BELs, both in their integrated line intensities and in their pro- 
files. In image A, we find that the profile of the C III] line is 
broadened during the brightening episode at HJD ~ 3900 days. 
The C III] line profile in image C seems also to be broadened 
at several epochs. Broadening of the BELs in image B is less 
obvious. 

Variations in the BEL intensities are detected mainly in im- 
ages A and B. Our measurements suggest that higher ioniza- 
tion BELs like C IV, C III], are more magnified than lower ion- 
ization lines like Mg II. This is compatible with reverberation 
mapping studies and a stratified structure of the BLR. There is 
marginal evidence that regions of different sizes are responsible 
for the Fe II+III emission. 

Finally, we estimate the differential extinction between 
pairs of quasar images due to dust in the lensing galaxy to be 
in the range 0.1-0.3 mag, with images C and D being the most 
reddened. This amount of differential extinction is too small to 
explain the large microlensing-magnification ratios involving 
images A and C. Long-term microlensing, beyond the duration 
of our observations, is therefore present in these images. 

The timescales of the microlensing variations in 
QSO 2237+0305 are such that each microlensing event 
lasts about one observing season (i.e. 8 months), with gaps of 
several months between events. This means that a relatively 
loose observing rate of 1 spectrum every 15 days is sufficient to 
sample the events well enough. In addition, the Einstein Cross 
is the lensed quasar with the fastest and sharpest microlensing 
events. It is therefore unique in the sense that only a few 
years of monitoring can truly constrain the quasar structure on 
parsec scales (Kochanek 2004). In addition, in the case of the 
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Einstein Cross, the very different behaviours of the BELs and 
the continuum with respect to microlensing offer considerable 
hope to reconstruct the two types of regions independently, 
using ray-shooting simulations. 

With two more years of data, we expect to map a total of up 
to half a dozen microlensing events in the four quasar images, 
hence providing a unique and useful data set for microlensing 
and quasar studies. 
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